Absolute determination of the 22 Na(p, 7) 23 Mg reaction rate in novae 
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Gamma-ray telescopes in orbit around the Earth are searching for evidence of the elusive ra- 
dionuclide 22 Na produced in novae. Previously published uncertainties in the dominant destructive 
reaction, 22 Na(p, 7) 23 Mg, indicated new measurements in the proton energy range of 150 to 300 
keV were needed to constrain predictions. We have measured the resonance strengths, energies, and 
branches directly and absolutely by using protons from the University of Washington accelerator 
with a specially designed beamline, which included beam rastering and cold vacuum protection of 
the 22 Na implanted targets. The targets, fabricated at TRIUMF-ISAC, displayed minimal degra- 
dation over a ~ 20 C bombardment as a result of protective layers. We avoided the need to know 
the absolute stopping power, and hence the target composition, by extracting resonance strengths 
from excitation functions integrated over proton energy. Our measurements revealed that resonance 
strengths for E v = 213, 288, 454, and 610 keV are stronger by factors of 2.4 to 3.2 than previously 
reported. Upper limits have been placed on proposed resonances at 198-, 209-, and 232-keV. These 
substantially reduce the uncertainty in the reaction rate. We have re-evaluated the 22 Na(p, 7) re- 
action rate, and our measurements indicate the resonance at 213 keV makes the most significant 
contribution to 22 Na destruction in novae. Hydrodynamic simulations including our rate indicate 
that the expected abundance of 22 Na ejecta from a classical nova is reduced by factors between 1.5 
and 2, depending on the mass of the white-dwarf star hosting the nova explosion. 

PACS numbers: 29.30.Kv, 29.38. Gj, 26.30. Ca, 27.30. +t 



I. MOTIVATION 

A classical nova is the consequence of thermonuclear 
runaway on the surface of a white-dwarf star that is ac- 
creting hydrogen-rich material from its partner in a bi- 
nary system. Such novae are ideal sites for the study of 
explosive nucleosynthesis because the observational [lj, 
theoretical @, Q , and nuclear-experimental 0, Q aspects 
of their study are each fairly advanced. In particular, 
due to the relatively low peak temperatures in nova out- 
bursts (0.1 < T < 0.4 GK), most of the nuclear reactions 
involved are not too far from the valley of beta stability 
to be studied in the laboratory, and the corresponding 
thermonuclear reaction rates are mostly based on exper- 
imental information 

An example of a gamma-ray emitter produced in novae 
is 26 Al (t 1/2 - 7.7 x 10 5 years, E 1 = 1.809 MeVF Indeed, 
this isotope has been observed in the Galaxy Q, but its 
long half life precludes the identification of its progenitor, 
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and novae are only expected to make a secondary contri- 
bution to its Galactic abundance 0, H| . Other gamma- 
ray emitters can provide more direct constraints on nova 
models [9]. An example is 22 Na (ti/2 = 2.603 years, E~ t 
— 1.275 MeV), which has not yet been observed in the 
Galaxy. Unlike 26 Al, the relatively short half life of 22 Na 
restricts it to be localized near its production site. Novae 
also could be the principal Galactic sites for the produc- 
tion of 22 Na, making 22 Na an excellent nova tracer. An 
observational upper limit of 3.7 x 10 -8 Mq was set on the 
22 Na mass in ONe nova ejecta with the COMPTEL tele- 
scope onboard the CGRO [lj3]. Currently, the maximum 
22 Na mass ejected using ONe nova models is an order 
of magnitude below this limit @, [H| and corresponds to 
a maximum detection distance of 500 parsecs using an 
observation time of 10 6 s with the spectrometer SPI on- 
board the currently-deployed INTEGRAL mission [12j . 
This suggests we are now on the verge of being able to 
detect this signal. In addition, Ne isotopic ratios in some 
meteoritic presolar graphite grains imply the in-situ de- 
cay of 22 Na produced by nucleosynthesis in novae fpjfTij 
and/or supernovae [Hf]. 

It is important to reduce uncertainties in the rates 
of key reactions that are expected to affect the pro- 
duction of 22 Na so that accurate comparisons can be 
made between observations and models 0]. For exam- 
ple, the production of 22 Na in novae depends strongly 
on the thermonuclear rate of the 22 Na(p, 7) 23 Mg reac- 
tion [1, 0, IHj], which consumes 22 Na. The thermonu- 
clear 22 Na(p, 7) reaction rate in novae is dominated by 
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narrow, isolated resonances with laboratory proton ener- 
gies E p < 300 keV. Consequently, the rate is dependent 
on the energies and strengths of these resonances, which 
have been investigated both indirectly and directly in 
the past. Indirect information on potential 22 Na(p, 7) 
resonances has been derived from measurements of the 
24 Mg(p,d) 0, 25 Mg(p,t) pjj, and 22 Na( 3 He,d) re- 
actions, and from the beta-delayed proton- and gamma- 
decays of 23 Al [13-113 ■ The nrst published attempt 
to measure the 22 Na(p, 7) reaction directly employed a 
chemically prepared, radioactive 22 Na target and pro- 
duced only upper limits on the resonance strengths [23j. 
A measurement contemporary to Ref. [H[ in the range 
E p > 290 keV by Seuthe et al. employed ion-implanted 
22 Na targets [24j, resulting in the first direct observation 
of resonances and the only absolute measurement of res- 
onance strengths. Later, Stegmuller et al. [H[ discovered 
a new resonance at 213 keV and determined its strength 
relative to the strengths from Ref. (24j. More recently, a 
new level in 23 Mg (E x = 7770 keV) has been discovered 
using the 12 C( 12 C,n7) [2(| reaction. This level corre- 
sponds to a 22 Na(p, 7) laboratory proton energy of 198 
keV, and the authors of Ref. [!|| proposed that this po- 
tential resonance could dominate the 22 Na(p, 7) reaction 
rate at nova temperatures. 

We have measured the energies, strengths, and 
branches of known resonances 124 , |25| and searched for 
proposed resonances [20Tj22l . |26| in the energy range 
Ep ~ 195 to 630 keV. The measurements were per- 
formed at the Center for Experimental Nuclear Physics 
and Astrophysics of the University of Washington with 
ion-implanted 22 Na targets prepared at TRIUMF-ISAC. 
Thanks to evaporated protective layers [27|, the targets 
exhibited little to no degradation over ~ 20 C of bom- 
bardment. Using mainly the strengths and energies ob- 
tained in this work to geth er with supplemental informa- 
tion from other work [24 l28j , we have re-evaluated the 
thermonuclear reaction rate of 22 Na(p, 7), and full hy- 
drodynamic simulations have been performed to estimate 
the effect of the new rate on the flux of 22 Na from no- 
vae. This article is a detailed presentation of our experi- 
ment, its results, and their implications, complementing 
our previous reports [H, Hfll • 



II. EXPERIMENTAL SETUP 

We measured 22 Na(p, 7) resonances directly by bom- 
barding implanted 22 Na targets with protons from a tan- 
dem Van de Graaff accelerator. High currents (~ 45 fj,A) 
at lab energies ranging from 150 to 700 keV were achieved 
with a terminal ion source. 



with areal density cINt / dA is given by 



prod 



dN T dN b 
~dA~dA 



dA, 



(1) 



where a is the cross section. Conventional methods em- 
ploy a small-diameter beam that impinges on a large-area 
target, where the target density is nearly uniform. How- 
ever, this technique can lead to target damage in cases 
where large beam currents are used, and there is a long 
history of differing results on resonance st rengths that 
have been attributed to target instabilities (3ll-l34|. We 
designed our experiment closer to the opposite limit, sim- 
ilar to Ref. 35J , where the beam was swept over an area 
larger than the full extent of the target with a rastering 
device. In the limit of uniform beam density over the 
target area, Eq. Q] becomes 
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This method requires knowledge of only the total number 
of target atoms and, thus, is not very sensitive to target 
non-uniformities. On the other hand, this method also 
requires a determination of the beam density. The yield 
is given by 



Y = (jN TPb , 



(3) 



where p b — ^j-/(Q/ e ) is a beam density normalized to 
the accumulated charge, Q. 

In addition, we determined the integrated yield of the 
excitation function over the beam energy, minimizing un- 
certainties associated with the energy loss in the target 
and beam energy distribution, which can be substantial 
in determinations using only the yield at a particular 
energy. The latter method, which was used in Ref. |24j . 
depends on knowing the energy loss in the target, the tar- 
get stochiometry and uniformity, and often assumes sta- 
ble target conditions, which are unlikely in experiments 
with currents of tens of microamps, such as ours. 

Beginning with Eq. [3] (see, for example, Ch. 4 of 
Ref. [3a]), the integrated yield for a finite-thickness target 
is given by 
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where J Y dE is the integral over the laboratory beam 
energy E with a range spanning the resonance, A is the 
reduced de Broglie wavelength in the center of mass, m 
and M are the projectile and target mass, respectively, 
and is the resonance strength. 



A. Strategy 

The number of reactions, A pro( j, produced by a beam of 
incident particles with areal density dN b j dA on a target 



B. Targets 

Our 22 Na targets were produced by ion implantation, 
which yields isotopically pure targets and avoids compli- 
cations with chemical fabrication. Each target was made 
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at TRIUMF-ISAC by implanting a 10-nA, 30-keV 22 Na+ 
ion beam into a rectangular OFHC copper substrate with 
dimensions 3 mm x 19 mm x 25 mm. The beam was 
rastered over a 5-mm collimator such that at the raster 
extreme only 5% of the beam remained on target, thereby 
creating a nearly uniform density. The setup included 
electron suppression with —300 V and a liquid nitrogen- 
cooled cold trap with a vacuum pressure in the range of 
6 x 10~ 8 to 2 x 10~ 7 torr. Charge integration was mon- 
itored throughout the implantation process, which took 
roughly 24 hours per target with a peak 22 Na current of 
- 15 nA. 

Initially, two test targets, #1 and #2, were implanted 
with activities of ~ 300 and 185 /xCi, respectively. As 
target degradation can be quite problematic, we carried 
out a program [27j to determine the ideal combination of 
implantation energy, substrate, and possible protective 
layer by bombarding 23 Na targets implanted under sim- 
ilar conditions. Using the conclusions of Ref. [27]], two 
additional ~ 300 ^Ci 22 Na targets, #3 and #4, were im- 
planted with the same parameters but included a 20-nm 
protective layer of chromium, deposited by vacuum evap- 
oration after implantation. A small rise in temperature 
of the target was observed during the evaporation; how- 
ever, a survey of the apparatus showed no residual activ- 
ity from diffusion of 22 Na out of the target. All 22 Na data 
presented were taken on the chromium-covered targets, 
with the exception of the 232-keV resonance measure- 
ment that employed target #1. 

To explore the transverse location of the implanted 
22 Na, the beta activity was scanned with a Geiger counter 
behind a 6-mm thick brass shield. A 3-mm diameter hole 
in the center allowed transmission of the beta particles. 
This measurement confirmed that the 22 Na was confined 
within a 5-mm diameter circle and determined the posi- 
tion of the activity relative to the center of the substrate. 
Although this method was not very sensitive to radially- 
dependent non-uniformity, it did verify axial symmetry. 
Thanks to this information and the extreme rastering of 
the 22 Na beam, we believe the targets were quite uniform, 
although our method does not require this. 



C. Chamber and Detectors 

Figs. [T] and [2] illustrate our chamber and detector sys- 
tem, respectively. The main features of the chamber in- 
cluded its dual liquid-nitrogen cooled cold shroud system, 
three collimators, and water-cooled target mount. The 
cold shroud isolated the radioactive target to prevent the 
contamination of the upstream beamline with 22 Na and 
also helped maintain a clean environment near the target, 
suppressing carbon buildup. During data collection, the 
pressure in the chamber was in the range of (1-2) xlO -7 
torr. The end of the downstream cold shroud surrounded 
our target substrate; however, because it was the farthest 
away from the liquid-nitrogen cold trap, it only reached 
a temperature of 125 K, whereas the upstream shroud 
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FIG. 1: (color online) Side view of chamber cross section. 
Copper braids connecting the upstream and downstream cold 
shrouds have been omitted for clarity. 



Beam 




FIG. 2: Top view of detector setup. The PMT shown is for 
the annular scintillator only; the planar scintillator PMT has 
been omitted for clarity. 



reached 88 K. 

The chamber had three sets of collimators. The first, 
collimator 1, was a water-cooled, sliding ladder between 
the cold shrouds with 4-, 7-, and 8- mm diameter collima- 
tors. The 8-mm collimator was used during 22 Na(p, 7) 
data acquisition. Also on this ladder were electron- 
suppressed 1-mm and 3-mm diameter collimators with 
downstream beam stops for tuning. Collimator 2, an 8- 
mm collimator, was 33 mm downstream of the ladder and 
was attached to the end of the downstream cold shroud. 
It was followed by a 10-mm diameter cleanup collimator, 
collimator 3, located 122 mm farther downstream. Each 
collimator was electrically isolated from the chamber to 
permit current monitoring. 
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The target substrate was bolted to a copper backing 
flange and was directly cooled with deionized-water via 
a thin pipe coupled to the flange. To minimize handling 
time in proximity to the radioactive targets, this assem- 
bly was then attached to a stainless steel coupler with 
a ISO flange on the chamber side. Directly upstream of 
this assembly was a 30-mm long electron suppressor bi- 
ased between —150 and —300 V. During data collection, 
the current on target was monitored, and the charge was 
integrated and recorded. 

Two sets of detector systems were positioned at ± 55° 
to the beam axis. A top view of the setup, including 
chamber, shielding, detectors, and Dewars, is shown in 
Fig. [2j Each system consisted of a high-purity 100% ger- 
manium (HPGe) crystal Canberra model GR10024 en- 
cased in cosmic-ray anticoincidence shielding. The detec- 
tor angle was chosen to minimize effects due to angular 
anisotropy, as the Legendre polynomial P2(cos8) is zero 
at ± 55°. The resolutions for each detector at 1.275 MeV 
were 4.4 keV and 7.4 keV (FWHM) with high rates ( 22 Na 
target present), and 2.2 keV and 3.0 keV (FWHM) with 
low rates (using residual 22 Na activity with 22 Na target 
removed). 

In addition, because of the target radioactivity, 26 mm 
of lead shielding was placed between the target and de- 
tector system to reduce the event rate in the detectors to 
a few tens of kHz. According to simulations, described 
in detail in Sec. HIE H the lead reduced the counting rate 
for the 511-keV gamma ray by a factor of 70, whereas 
the photopeak from 1275-keV gamma rays was reduced 
by a factor of 5. The suppression ranged from factors 
of ~ 3.5 to 4.5 for gamma rays with energies above ~ 
5 MeV. Above 4 MeV, cosmic rays caused a continuum 
background. In order to remove these unwanted signals, 
25-mm thick annular and planar plastic scintillators en- 
cased in lead, as shown in Fig. [51 were used in anticoinci- 
dence with the germanium detectors and are discussed in 
detail in the next section. The reduction of the 22 Na(p, 7) 
signal by this veto was negligible. 



D. Data Acquisition 

In addition to the signals from the high-purity Ge de- 
tectors, the electronics also processed PMT signals from 
the scintillators. In order to reduce the rate seen by the 
detectors due to target radioactivity, two techniques were 
used: 1) lead shielding was installed as described above, 
and 2) a high threshold was set (just below the strong 
1275-keV 22 Na line, so that the target activity could be 
monitored in-situ). 

The Ge-signal amplifiers (ORTEC 672) were operated 
in a pile-up-rejection mode (which typically rejected ~ 
40% and introduced a dead time of 27 [is per pulse). 
Signals from the two sets of detectors were converted 
into digital signals by ORTEC AD413A ADCs with fast 
FERAbus readout, which helped to reduce dead time. 
A buffer module was used to minimize the communica- 
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FIG. 3: Sample background spectrum from one Ge detec- 
tor, measured for ~ 160 hours. Inset shows the spectrum 
expanded around the 1275-keV gamma ray. 

tion with the computer via a CAM AC interface. JAM, 
a JAVA-based data acquisition and analysis package for 
nuclear physics [37j, was used to process the data. All 
the NIM and CAMAC electronics modules were located 
in a temperature-controlled rack to minimize instabili- 
ties. The raw rate in each detector was below 30 kHz, 
and the trigger rate was ~ 4 kHz. A sample background 
spectrum from one Ge detector is shown in Fig. [3] 

Signals from the active anticoincidence shields and 
from the Ge detectors were used as a stop and start, 
respectively, in a Time-to- Amplitude- Converter (TAC). 
If these two signals occurred within a set timing window, 
the resulting Ge signal was discarded. 

In order to determine the bounds of the TAC spec- 
trum for data processing, the TAC signal corresponding 
to germanium detection energies between 4 and 6 MeV 
was extracted. The timing gate was set on the prompt 
peak, which had a long tail to its left. In this energy re- 
gion, this anticoincidence system rejected ~ 80% of the 
cosmic-ray background signal. The TAC signal for most 
high-energy cosmic rays from ~ 6 to 11 MeV also falls 
within the set window. The scintillator threshold was set 
well above 511 keV to avoid vetoes by annihilation radi- 
ation. Self veto is possible with cascade gamma decays if 
one gamma ray is registered in the Ge and the other in 
the anticoincidence shield; this was examined by compar- 
ing the yields from raw singles spectra with those from 
the vetoed spectra, which agreed to better than 99%. 



E. Detector Efficiency 

Detector photopeak efficiencies were obtained by com- 
bining Penelope (3|1 139[ simulations with direct mea- 
surements of gamma rays from 60 Co and 24 Na sources 
and from 27 Al(p, 7) resonance measurements. 

The efficiency at 1332 keV was determined using a 
31.51-nCi 60 Co source. Isotope Products produced 
the source and measured its decay intensity with an un- 
certainty of 1.7% (99% C.L.). The geometry for the 



5 



Penelope simulations was adjusted slightly to match 
the efficiency at this energy, as described in the following 
subsection. The ratio of efficiencies from 1369 to 2754 
keV was measured using a 24 Na source (t 1 / 2 = 15 hrs) 
fabricated at the University of Washington. Since 1369 is 
very close to 1332 keV, we equated the efficiency at 1369 
keV to the Penelope value, and then from this ratio we 
obtained the efficiency at 2754 keV. 

To extend the efficiency determination to higher ener- 
gies, we measured the 27 Al(p, 7) reaction using a thick 
aluminum target and used the relative intensities of 
well-known gamma-ray branches from E p = 633 and 
992 [U HH. For the latter resonance, our coin target 
was also used. We subtracted the below- resonance yield 
from the above-resonance yield to extract the net yield. 
At E p — 992 keV, the gamma rays of interest are at 
1779, 4742, and 10762 keV. Using the simulation to de- 
termine the efficiency at 1779 keV, our measurements, 
combined with the known branches, gave the efficiencies 
at the two higher energies. At E p — 633 keV we measured 
the gamma rays at 7575 keV and 10451 keV. As the sim- 
ulation matched the value we had obtained at 10762 keV, 
we used the simulation for the 10451 keV value and the 
known branch to determine the efficiency at 7575 keV. 
The agreement between measurement and simulation is 
discussed in the following subsection. 



1. Penelope Simulations 

The geometry of our apparatus was modeled in the 
detailed Monte Carlo code Penelope. The simulated 
germanium detector included the germanium crystal, 
cold finger, and carbon window, with all dimensions 
taken from the nominal specifications provided by Can- 
berra [43]. In addition, we included the 26-mm lead and 
25-mm planar plastic scintillator in front of the detector. 
Although the annular plastic scintillator was modeled, 
the annular lead was not, as it was not in the line of 
sight of the target. The sodium source was a uniform 
5-mm diameter disk centered on the copper substrate. 
The copper backing mount was included, as was the alu- 
minum plate supporting the water cooling system. The 
water and its copper pipes inside the target mount were 
modeled, but the pipes that extended up and out from 
the mount were not, as they were thin and mostly out 
of the line of sight. All components of the target mount 
were aligned with the beam, whereas the detector was at 
an angle of 55°. 

The gamma rays were projected from their source uni- 
formly in a 80° opening angle, which covered all modeled 
components, and absolute efficiencies were corrected for 
the solid angle. Each simulated energy included an ini- 
tial number of gamma rays such that the photopeak pre- 
cision was less than 0.1%. At E 1 =1332 keV, with the 
source spread out over an area equal to that of the 1- 
mm diameter 60 Co source, the simulation initially gave 
results 2.5% higher than the measurement. Therefore, 



7x10" 4 

6 - 

a 

.2 4- 

o 

5-< 7 _ 

2 H 

1- 
1.05 
■$ 1.00 
£ 0.95 
0.90 



- simulation 

* 6 °Co 
» 24 Na 

> 27 A1( P ,7): 
E p =992 keV 

• 27 Al(p, 7 ): 
E =633 keV 

■ 27 A1(P,7): 
E =992 keV 



6 8 

Energy (MeV) 



10 



FIG. 4: (color online) Photopeak efficiency. Top panel is 
the efficiency from Monte Carlo Penelope simulations. Bot- 
tom panel shows the ratio of efficiencies from measurement to 
simulation, and arrows indicate the gamma-ray energy used 
in the relative efficiency analysis. 



the front face of the crystal was moved back from the 
target by 1.7 mm, in order to make the simulation re- 
produce the measurement exactly. Results for the detec- 
tor photopeak efficiency are shown in the upper panel of 
Fig. [4] The bottom panel of Fig. [4] shows the ratio of 
the measured efficiencies to the simulations. For sources 
other than 22 Na, source distribution and substrate mate- 
rial were changed in the simulation to match those used 
in the measurement. 

For the 213- and 610-keV resonances in 22 Na(p,7), 
yields from first-escape peaks were added to the pho- 
topeak yield in order to improve statistics for branches 
with E 1 = 7333 and 8162 keV, respectively. Compari- 
son of 27 Al(p, 7) data to simulation at 7.5 MeV indicates 
agreement to within 2% and is covered by the systematic 
error detailed in the next section. 



2. Systematic Errors for the Efficiency 

In order to extract systematic errors for our efficien- 
cies, we compared the quality of the fit of our data from 
the 24 Na source and the 992- and 633-keV resonances of 
27 Al(p, 7) to our simulations, as shown in Fig. 2] The pre- 
cisely measured ratio for the two 24 Na gamma rays yields 
a value for the efficiency at 2754 keV which is in agree- 
ment with that given by simulation. The points obtained 
from 27 Al(p, 7) resonances have statistical uncertainties 
between 4.7 and 5.2%, and they agree well with the simu- 
lation. Therefore we ascribe a 5% systematic uncertainty 
to the efficiency determination for isotropic emission of 
gamma rays. 

Because the detectors are centered at ± 55° in the 
laboratory, zeros of /-^(cosfl), the effect of a P 2 (cos6) 
term in the angular distribution can only arise from the 
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angular dependence of the efficiency across the detec- 
tor and from center of mass to laboratory transforma- 
tion. A P4(cos6>) term has no such restriction. Assum- 
ing the angular distribution to be of the form = 
1 + a 2 P2(cos#) + a4P4(cos#), we used the Penelope 
simulation to determine the effect of non-zero values of 
a-i and 04. A value of ai as large as 1 only caused a 
2.6 ± 0.4% change in the efficiency. Published data for 
23 Na(p, 7) resonances [44[ show typical 04 values of about 
0.005 and a maximum value of 0.05. This maximum value 
would cause a 2.0 ± 0.4 % change in the efficiency. There- 
fore, we assigned an additional systematic error of 3%, 
to include the possible effects of the angular distribution. 
Our overall systematic error in the efficiency is ± 6%. 

F. Beam Properties 

1 . Rastering 

The beam was rastered using a magnetic steerer lo- 
cated f m upstream of the target. A rectangular pattern 
with 19- and 43-Hz horizontal and vertical frequencies 
was used, and signals proportional to the magnetic field 
were produced by integrating the voltage signals from a 
pickup coil located in the raster magnet. These read- 
out values represent the center of the beam spot. For 
each data set collected, a 2-dimensional histogram of 
this signal in both horizontal and vertical directions was 
recorded. It was also possible to set a gate on the energy 
spectrum of each detector and sort out the corresponding 
raster fields. 

Fig.[5]shows this two-dimensional raster plot, obtained 
with a 5-mm diameter, 1.5-mm thick 27 Al disc embedded 
in the center of a copper backing. This "coin" target had 
the same OFHC copper substrate and diameter as our 
22 Na targets. Shown are the counts detected for each 
value of the rastering field in the two dimensions. Dur- 
ing all data-taking, the raster signals were monitored in 
order to diagnose problems with the target or other issues 
that may have arisen. For each measured resonance, the 
amplitude for the raster was scaled as the square root of 
the proton energy, so the rastered area of the beam on 
target would be constant. 

2. Energy and Density 

The beam energy was determined using a 90° analyzing 
magnet with an NMR field monitor and calibrated using 
resonances in 27 Al(p,7) at E p = 326.6, 405.5, 504.9, and 
506.4 keV. Given the quality of the fit, we assign a ± 0.5 
keV uncertainty to our knowledge of the beam energy. 

We also conducted measurements on 27 A1 targets in 
order to extract the normalized beam density, averaged 
over the 5-mm diameter area. Systematic errors were 
then explored with simulation and visual inspection of 
beam-related target coloration. Beam non-uniformity 
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FIG. 5: (color online) 2-dimensional raster plot on the "coin" 
target: a 5-mm diameter 27 A1 disc embedded in copper. Bx 
and By are proportional to the magnetic field of the raster. 
Shown is a histogram of the raster signal during a measure- 
ment of the 406-keV 27 Al(p, 7) resonance, with a wide gate 
ranging from above the 7358-keV gamma ray to below its 
first-escape peak. 



was then taken into account, and its effect combined with 
the possibility of a non-uniform target was investigated. 

We determined the normalized beam density, pb, aver- 
aged over the target, by comparing the thick-target yield 
from our "coin" target, Y c , which had the same areal ex- 
tent as the 22 Na targets, to the yield from a solid 27 A1 
target, Y s . From the ratio of these yields, one can extract 
p b via: 



7^ = PbA c , (5) 

where A c is the area of the coin. The yields were mea- 
sured at different times with different beam tunes using 
the resonances at E p = 406 and 992 keV, which yielded 
results for p b of 2.60 ± 0.09 cm" 2 and 2.55 ± 0.10 cm' 2 , 
respectively, with statistical errors only. The weighted 
average of these two measurements, 2.58 cm~ 2 , was cho- 
sen for pb- 

To study the distribution of beam across the target 
in order to quantify systematic errors, we carried out a 
number of measurements using a large raster with 1.8 x 
the standard amplitude, in addition to the standard am- 
plitude, on the main 22 Na targets and the solid 27 A1 
and coin targets. With the raster on and off, transmis- 
sion measurements through various collimators, shown 
in Fig. [TJ were performed as well. These measurements, 
along with the relative yields of large to standard raster, 
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FIG. 6: The relevant part of the decay scheme of 23 Mg. En- 
ergies are in keV. Not all 23 Mg levels are shown. 



were used to constrain a Monte Carlo simulation, de- 
scribed below, that was used to investigate potential 
beam densities on the target. It also was used to test 
the effects of possible beam drift, misalignment of the 
beam and target, and target non- uniformities. 

The simulation modeled transport of the beam through 
the final components of the beamline and chamber, in- 
cluding the final quadrupole, the rastering unit, and the 
three sets of collimators. Variable parameters included 
beam width and offset, possible collimator offset, raster 
amplitude, and beam distribution at the quadrupole. A 
normalized beam density could not be uniquely deter- 
mined by this method alone, but the densities within an 
acceptable phase space, defined by reasonable agreement 
for transmissions and large/standard rastering yield ra- 
tios, ranged from 2.33 to 2.83 cm~ 2 . Even for the extreme 
case where the quadrupole aperture is filled uniformly, 
the beam density was found to vary by less than 15%. 
We adopted ± 0.25 cm" 2 as the systematic uncertainty 
in the normalized beam density. 



III. DATA AND ANALYSIS 

Measurements were taken on previously known 
22 Na(p, 7) resonances, which we find at E p = 213, 288, 
454, and 610 keV. We also explored the proposed res- 
onances at 198, 209, and 232 keV. The relevant energy 
level diagram for 23 Mg is shown in Fig. [SJ 

For all measurements, an energy scan was performed 
across a ~ 25 keV range about the nominal resonance 
energy. To subtract background, which was comprised 
mostly of cosmic rays and Compton events, we assumed 
it to have a localized linear dependence on gamma-ray 
energy and fit the background to windows in the spec- 



trum above and below the line of interest. This method 
was especially important for resonances with charac- 
teristic gamma rays below 6 MeV. Here contributions 
from 6129-keV incident photons due to the contaminat- 
ing 19 F(p, a7) 16 reaction were significant. Affected res- 
onances included E p = 454, 288 keV, and one branch 
from E p = 610 keV. 

In order to extract the yield at each laboratory proton 
energy, E p , an energy window was set on the gamma ray 
of interest in the vetoed singles spectrum for each of the 
two germanium detectors. This window was ~ 25 keV 
for one detector and ~ 40 keV for the other. For each 
detector, the sum of the background subtracted counts in 
the window is Ni = Np^rjiLi, where rji is the efficiency, 
and Li is the live time for detector i. Then the yield is 



Y 



jVprod = A1/L1 + N 2 /L 2 
Q/e 



r)Q/e 



(6) 



where r\ = 771 + 772. The effects of angular distributions 
have been addressed in Sec. Ill E 21 In order to deter- 
mine the live time, a signal from a pulser unit was fed 
into the "test" port of each Ge preamplifier, creating an 
additional signal in the corresponding amplitude spec- 
tra. This signal was sorted into its own spectrum by a 
logic signal to the data acquisition. A window of com- 
parable width to the energy window for the yield was 
placed on the prompt pulser signal and compared to the 
scaled number of pulser pulses. This ratio gave the live 
time, which ranged from 35 to 45% for the radioactive 
targets and was above ~ 90% for all other targets. The 
live-time correction was substantial, but it was not beam- 
related; instead it resulted from a constant rate due to 
the radioactivity. To test the accuracy of the live-time 
correction, a thick 27 A1 target was irradiated with pro- 
tons with and without a 22 Na source nearby. Although 
the presence of the radioactivity decreased the live time 
from 97% to 50%, the ratio of the live-time corrected 
yields was 0.99 ± 0.02. 

In order to test the sensitivity of our results to the in- 
puts for the gamma-ray background subtraction and its 
linearity, resonances at E p — 610 keV (with E 7 = 8162 
keV) and 454 keV (with E 1 = 5300 keV) were inspected, 
and the choice of window for both the peak and the 
background on each side was varied to reasonable limits, 
such as widening, shortening (no window smaller than 10 
keV), and shifting. The resonance strengths changed by 
less than 1%, indicating that systematic errors associated 
with background subtraction are negligible. 

After the yields for each excitation function were de- 
termined, the areas under the excitation curve were es- 
timated by using the trapezoidal method. This value 
for J Yi dE was then used in Eq. 0] along with values 
determined for all other parameters, to extract the par- 
tial resonance strength, uryj, for each branch i. The to- 
tal strength is simply equal to the sum of the partial 
strengths for all branches. 
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FIG. 7: Excitation functions for (a) E p = 454 keV and (b) 
E p = 610 keV. E p is the lab proton energy. Each plot shows 
the excitation function at the beginning of target bombard- 
ment and at the end, after ~ 20 C had been deposited. An 
intermediate curve after 10.5 C is also shown in (b). (a) is 
gated on E 7 = 5300 and 5963 keV, and (b) is gated on the 
photopeak and single-escape peak of E 7 = 8162 keV. 



FIG. 8: Summed gamma-ray spectra for (a) Ep = 454 keV and 
(b) Ep = 610 keV. The inset illustrates the beam-background 
subtraction method for E 7 = 5300 keV, including gamma-ray 
gate (dashed), analytically calculated background line (solid), 
and subtracted region (hatched). SE and DE indicate single- 
and double-escape peaks, respectively. 



A. Yields 

1. Absolute Yields: E p = 610, 454, 288, 213 keV 

Fig. [7] shows the data taken on the two strongest reso- 
nances at E p = 454 and 610 keV. These resonances were 
revisited after various amounts of accumulated charge to 
monitor possible target degradation, discussed in detail 
in Sec. 1111 Bl Fig. [5] shows the corresponding gamma-ray 
spectra summed over all runs, including an inset illus- 
trating the background subtraction method. All data for 
resonances at E p = 213, 288, and 610 keV were taken on 
target #4, and the 454-keV resonance was measured on 
both targets #3 and #4. 

Fig. IH1 illustrates the summed raster plots for target #3 
with E p = 454 keV. The concentration of counts from 
22 Na(p, 7) are well centered, while a few counts spread 
through the plot are consistent with yield from 19 F con- 
tamination. 

Fig. [TCil shows the data taken for 288- and 213-keV res- 
onances, and Fig. [TT] illustrates the summed gamma-ray 
spectra for each resonance, respectively. Characteristic 
gamma rays for each are clearly distinguishable above 
background. For E p = 213 keV, the integrated yield was 



determined with the same trapezoidal method as other 
resonances, plus a small correction because the highest 
energy point did not reach zero yield. The details of de- 
termining this contribution are discussed in Sec. 1111 A 31 
after a prerequisite analysis method is outlined in the 
following section. 



2. Upper Limits for Yield: E v = 198, 232, 209 keV 

Data in the region of the proposed resonances at 
E p = 198 and 232 keV are shown in Fig. [H In the 
summed gamma-ray spectra shown in Fig. 1131 no dis- 
cernible gamma-ray yields can be detected above back- 
ground. All data for E p = 198 keV were taken on the 
chromium-covered target #3. Data for E p = 232 keV 
were taken on one of the bare test targets, which had 
previously been exposed to an integrated charge of 13 C. 

For these resonances, the shape of the excitation func- 
tion used to determine the area was adopted from cither 
the resonance at E p = 454 or 610 keV, depending on the 
target. Because the resonance shape was dominated by 
the implantation distribution, this shape was normalized, 
shifted, and stretched so that it could be fit to the data 
of the resonance in question. The stretch factor was fixed 
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FIG. 9: (color online) Raster plot for the sum of 22 Na(p, 7) 
for E p — 454 keV, target #3, with raster amplitudes 1.8x 
larger than our standard amplitudes. The plot is gated on 
E-y — 5300 keV. Outside the target area there are a few bins 
with only 1 count, which is consistent with the known level 
of 19 F contamination from the Compton continuum from E 1 
= 6129 keV. 



and set equal to the ratio of stopping powers in copper for 
the two energies, whereas the energy shift and the nor- 
malization factor were allowed to vary. The central value 
of the shift was given by the differences in resonance en- 
ergies, and the range of the shift was given to fully span 
the data points. If a data point for the low-energy res- 
onance fell between points of the normalized curve, the 
corresponding reference point was determined by a linear 
interpolation. For each pair of shift and normalization, 
the value of the \ 2 between each low-energy resonance 
and the normalized reference resonance was calculated. 
The modified reference excitation function corresponding 
to the minimum value of \ 2 is shown for each resonance 
in Fig. [12] (open circles) . 

The array of probabilities, P(aryi, E{), was taken to be 
proportional to exp(— x 2 {^li, Ei)/2), where x^i^li^i) 
is the x 2 between the model, assuming particular val- 
ues for ujji and Ei, and the data. Because we are 
mainly interested in constraining the value of ory, we pro- 
jected the two-dimensional arrays onto the W7-axis (i.e. 
P(w7i) = YjE P{^H^ E i))i shown in Fig.QH The upper 
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FIG. 10: Excitation functions for (a) E p = 288 keV and (b) 
E p is the lab proton energy, (a) is gated 



Ep = 213 keV. 



on E 7 = 5141 keV, and (b) is gated on the photopeak and 
single-escape peak of E 7 = 7333 keV. 
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FIG. 11: Summed gamma-ray spectra for (a) E p — 288 keV 
and (b) E p = 213 keV. SE and DE indicate single- and double- 
escape peaks, respectively. 
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(color online) Solid circles are excitation functions 



E p — 



232 keV and (b) E p 



proton energy, (a) is gated on 



= 198 keV 
5055 keV, and (b) is gated 



E p is the lab 



on the sum of E 7 = 7801 and 7350 keV. Open circles are the 
normalized, stretched, and shifted E p = 454 keV excitation 
function fit, described in the text. This reference excitation 
function is not as smooth in (a) as in (b) due to degradation 
of the bare target. 



FIG. 13: Summed gamma-ray spectra for (a) E v = 232 keV 
and (b) E v = 198 keV. The target activity was 2.3 x higher in 
(b) than in (a) and is responsible for the increased background 
at low energies. Arrows indicate the energy where one would 
expect to see gamma rays from the relevant transition. 



limits on W7 were extracted with a particular confidence 
level, C.L., using the likelihoods: 

CL. = ^ . (7) 

The sum in the denominator was cut off at a maximum 
value of twy such that the sum changed by less than 1%. 
Results are given in Sec. HVl 

For the possible resonance at 198 keV, a finite value 
for the strength was also calculated. Instead of sum- 
ming from zero and extending upward in the numera- 
tor of Eq. [71 the pair of uiji values with equal values 
of P(w7i) were determined such that the sum between 
them, properly normalized, gave the desired confidence 
level. Because of its small branch, data for the third pos- 
sible gamma-ray for E p = 232 keV at E 1 = 5749 keV was 
not added to our yield; however, branches from the two 
other gamma rays were used to adjust the total resonance 
strength. 



For the proposed resonance at E„ = 209 keV, we also 
applied this method to a hybrid data set comprised of 
measurements around the resonances at E p — 198 to 213 
keV for the gamma ray at 5067 keV, assuming the branch 
given by Jenkins et al. [26| and using the excitation func- 
tion from the dominant branch of E p = 213 keV and its 
first-escape peak as the reference curve. Because the data 
from the resonance at 198 keV were from a different tar- 
get, its yields were scaled by the ratio of measured target 
activities. The shift in energy was allowed to vary from 
zero to 25 keV, and the fit yielding the minimum x 2 wa s 
found at the position of the 198-keV data points, possi- 
bly due to the fact that the gamma rays have overlapping 
energy windows. In other words, we did not observe a 
separate resonance at E p — 209 keV. An upper limit for 
this resonance, which is presented in Sec. IIV1 was ex- 
tracted by restricting the energy shift to be equal to the 
difference in resonance energies, spanning the r ang e of 
± 2<7 around the value claimed by Jenkins et al [26j. 

This analysis technique of normalizing and shifting a 
reference resonance to obtain strengths of others was val- 
idated by applying it to the 288-keV resonance, for which 
the ratio of the strength calculated from this method to 
the direct method was 0.95 ± 0.12. 
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FIG. 14: Projection of exp(— x 2 (c^7i, Ei)/2) onto the W7-axis 
for (a) E p = 232 keV and (b) E p = 198 keV. The curves 
represent the probability density functions, used in Sec. [V] 



3. Corrected Area for E p = 213 keV 

In order to estimate the full area of the E p = 213 keV 
excitation function, a reference resonance measurement 
at E p = 454 and three at 610 keV were utilized in the 
same manner outlined above. Each of the four curves 
were fit to the data, and each yielded a data point beyond 
the fixed E p — 213 keV excitation function that did reach 
zero. The last trapezoid area was calculated for each, and 
the average was added to the area from the direct data, 
equaling 10 ± 5% of the total area. The uncertainty in 
the additional area was set to be the standard deviation 
among the four fits. 



B. Total Number of Target Atoms 

We determined the initial number of target atoms from 
the 1275-keV gamma r ays emitted in the decay of 22 Na 
(ti/2 = 2.6027(10) yrs [45]) and assigned a 2.6% uncer- 
tainty. This uncertainty combines the 1.7% uncertainty 
in the 60 Co calibration source at E 1 = 1332 keV with 
an additional 2% due to the accuracy of the ~ 7% back- 
ground subtraction in this region of high detector rate. 

However, measuring the activity in-situ was not suf- 
ficient to determine the total number of target atoms 
throughout the measurement. During target bombard- 
ment, some fraction of 22 Na was sputtered out of the il- 




Charge Q (C) 

FIG. 15: (color online) Target degradation as a function of 
accumulated charge for main targets #3 and #4. Nt(Q) is 
the total number of atoms present in the substrate after an 
irradiation of charge, Q. The data points illustrate the ratio 
of the area of each reference excitation function to its initial 
area at the start of target bombardment. The lines represent 
a possible linear decrease of the number of atoms deduced 
from two residual activity measurements at the beginning and 
end of bombardment. Also shown is the timeline for each 
resonance measurement. The two sections marked "454T" 
indicate variable raster amplitude tests, and although these 
data are included in the branching ratio determination, they 
are not included in determining resonance strengths. 



luminated area of the substrate, yet it remained nearby, 
maintaining an approximately constant activity through- 
out the duration of the resonance measurements. Thus, 
in addition to determining the total number of initial 
atoms, monitoring possible target degradation through- 
out bombardment was particularly important. 

In order to deduce the amount of degradation, two 
complementary methods were utilized. One method 
was to revisit a strong reference resonance periodically 
throughout the bombardment cycle. We define A(Q) as 
the integral of the reference excitation function, J YdE, 
after an amount of charge, Q has been deposited. A(Q) 
is directly proportional to the number of target atoms, 
Nt(Q), as shown in Eq. 2) The ratio of the integrals of 
the excitation function before and after bombardment, 
A{Q)/A{Q - 0), is therefore equal to N T (Q)/N T (Q ~ 
0). The second method, which will be described in detail 
later in this section, was to measure the residual 22 Na in 
the chamber before and after target bombardment and 
use this information to infer the number of sputtered tar- 
get atoms. The results of each method are illustrated in 
Fig. 1151 along with a timeline for each resonance mea- 
surement as a function of accumulated charge. 

Target #3 accumulated 20.7 C in 186.0 hours, and tar- 
get #4 accumulated 20.1 C of charge in 138.3 hours. Be- 
cause we covered these targets with 20 nm of chromium, 
they remained fairly stable. Shown in Fig. [7] (a) are the 
first and last resonance scans at E p = 454 keV for tar- 
get #3. Throughout the 20.7 C of bombardment, the 



resonance at E v 



454 keV was revisited four times. 



A{Q ~ 20 C)/A(Q ~ 0), shown in Fig. [Tj3 is 1.07 ± 
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0.12, consistent with no target loss. For target #4 on 
which all other non-zero strength resonances were mea- 
sured, the monitoring resonance was E p = 610 keV, and 
multiple scans of its excitation function are shown in 
Fig. [7] (b). At the end of bombardment, A(Q ~ 20 
C)/A(C - 0) = 0.94± 0.09, which is consistent with 
no target loss. 

To use the residual activity method, measurements of 
the 1275-keV rate were taken before target installation 
(J?q), after target installation and before bombardment 
after bombardment (R T ), and after target removal 
(Rq). The quantity (R!^ — R l Q ) — (R^. — R( j ) then is propor- 
tional to the amount sputtered from the target. This final 
value was used to estimate target degradation through- 
out the bombardment, assuming linear loss, and is also 
shown in F ig. | 15l for each target. As we learned with our 
23 Na tests [27|, this loss is in fact not linear but usually 
begins to occur after a significant amount of charge has 
been deposited, removing the protective layer and sput- 
tering away some of the substrate. We nevertheless used 
a hypothesis of linear degradation of the target for one 
extreme of Nt and an amount consistent with no loss for 
the other. 

The 198-, 213-, and 232-keV resonance measurements 
were taken over an extended period of time and charge, 
whereas all others were measured with a few Coulombs 
of integrated beam current and did not experience pos- 
sible prolonged degradation. The 198- and 213-keV data 
were taken over ~ 15 and 10 C, respectively, with short 
interruptions to measure the reference resonance excita- 
tion functions. At the halfway point for each, the linear- 
decrease hypothesis indicated a 4-5% loss, although exci- 
tation function areas are consistent with no loss at that 
point. Therefore, we choose no loss with errors that 
span the values from each method, Nt(Q)/Nt(Q ~ 0) = 
1.00^0 05- Combining this with the systematic uncer- 
tainty in the initial number of atoms, we have an overall 
systematic error of +2.6% and —5.6% in the total number 
of atoms for the 198- and 213-keV resonance strengths. 
For the 454- and 610-keV resonances, which were each 
measured at the beginning of target bombardment, only 
an overall systematic error of ±2.6% was needed. 

The 288-keV resonance was a special case, as its data 
were not from an extended measurement but were taken 
after 18 C of irradiation. Directly following the measure- 
ment of this resonance, we performed the final scan of the 
610-keV reference resonance, which allowed insight into 
how many atoms remained. Therefore, the total number 
of atoms present for the 288-keV measurement was taken 
to be the average between linear target loss and loss in- 
dicated from the depleted area of the 610-keV resonance 
curve, N T {Q)/N T (Q - 0) = 0.88±g;J§. Because target 
loss is not actually linear and could have happened after 
the 288-keV resonance measurement, the uncertainties 
span the range between no loss and the value given by 
linear loss. 

The 232-keV resonance data were taken with a test tar- 
get with no protective layer after 13 C had already been 



TABLE I: Summary of systematic errors for each resonance 
strength. Note, the total error for the resonances at E p = 198 
and 232 keV is dominated by statistical errors. 
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bombarded, and target loss was appreciable. At the end 
of the ~ 20 C irradiation, measurements of the resid- 
ual activity indicated a 68% loss. As explained above, 
loss is not linear and occurs quite rapidly at the end of 
the cycle. Because of this fact, we have chosen to take 
the loss only from the difference in monitoring the reso- 
nance areas, which were taken directly before and after 
the measurement and gave A(Q)/A(Q ~ 0) = 0.59. This 
value was used to adjust the upper limit on the strength, 
and we assigned a systematic uncertainty of ± 40% to 
span a wide range approximately down to the value of 
linear loss. Regardless, the upper limit on this resonance 
strength is still dominated by the statistical uncertainty. 



C. Systematic Error Summary 

The systematic error budget for resonance strengths 
is shown in Table HI In summary, we find an overall 
systematic error of —11.7% and +12.7% for the extended 
measurements at E p — 198 and 213 keV, resulting from 
combining uncertainties of ±6% in the efficiency, ±10% 
in the normalized beam density, and —2.6% and +5.6% 
from the number of target atoms. For the resonances at 
E p = 454 and 610 keV, the overall systematic uncertainty 
is ±11.7%, and for the resonance at E p = 288 keV it is 
— 18.1% and +13.4%, differences all due to the differing 
number of target atoms. The 232-keV resonance has the 
largest overall uncertainty, ±42%. 



D. Resonance Energies 

We were able to obtain resonance energies with two 
separate techniques. First, we extracted the resonance 
energy from the observed gamma-ray energy, along with 
the excitation energy of the daughter level [45j and Q- 
value (7580.53 ± 0.79 keV, using the newly measured 
masses of 23 Mg HI] and 22 Na [47|). From a thick-target 
27 Al(p, 7) measurement at E p — 406 keV, the spectrum 
was calibrated using transitions energies 5088.05 and 
7357.84 keV and each first-escape peak from the cor- 
responding gamma rays. Energies were corrected for 
Doppler shift and recoil, which ranged in magnitude from 
~ 4 to 7 keV and ~ 0.6 to 1.5 keV, respectively. Because 
the detector gain depended slightly on rate, a correction 
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ranging in magnitude from 2.4 ± 0.7 to 3.8 ± 1.1 keV was 
also applied. This was determined by measuring the shift 
in the calibration gamma-ray energies with and without 
22 Na sources nearby. 

Second, we found the resonance energy from the proton 
energy, via the excitation function. The energy at which 
the yield reached half its maximum was determined, and 
the losses in the 20-nm chromium layer and 4 nm of cop- 
per were subtracted. This copper depth is the depth at 
which the 22 Na distribution reached half of its maximum 
value. According to simulations using TRIM (48|, the 
total subtracted losses were ~ 3 to 5 keV, and a 20% 
uncertainty in the stopping power was assumed. For the 
213- and 288-keV resonances, an additional adjustment 
to the resonance energy was added to account for the 
slight transformation of the excitation functions due to 
sputtering, as shown in Fig. 0(b). From repeated scans 
of the 610-keV resonance, the energy at half of the max- 
imum yield changed by 1.2 ± 0.9 keV after 11 C, in the 
middle of the 213-keV resonance measurement, and that 
shift remained constant after 19 C, directly after the 288- 
keV measurement. Those resonance energies were ad- 
justed by that amount. Results are given in Sec. IIVI 
Similarly, excitation energies, E x , and gamma-ray ener- 
gies, Ej, each were found independently by using the 
weighted average of the respective value extracted from 
the excitation function and the respective value extracted 
from the gamma-ray spectra. 
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FIG. 16: Excitation function for 23 Na(p, 7) resonance at E p 
= 512 keV gated on — 10810 keV. E p is the lab proton 
energy. 

E 1 = 8162 keV, the largest branch in the de-excitation. 
The value of this normalization was equal to the ratio 
of the magnitude of the branches, adjusted for efficiency 
differences, and was extracted by minimizing the x 2 of 
the hybrid curve and the on-resonance curve. Because 
of the changing shape of the off-resonance curve above 
and below the resonance, only an upper limit could be 
obtained. 



E. Branches 

Strong branches were determined from the spectra 
summed over all runs within a particular resonance. For 
resonances at E p — 454 and 610 keV, which were used 
as reference resonances to monitor degradation and for 
other target tests, the total amount of data was signifi- 
cantly larger than for a single resonance scan. 

Due to the very low statistics for weaker possible 
branches, an additional restriction was placed on the 
analysis. Similar to our analysis for weak resonances 
where a reference resonance excitation function was 
shifted and normalized to fit the data, the excitation 
function for the strongest branch (and its first-escape 
peak for branches from E p = 610 keV) was normalized 
to match the weaker branches' excitation functions, such 
that the y 2 was minimized. Then the branches were ex- 
tracted using Eq. [71 

Only an upper limit for the branch at E~ = 6112 keV 
from the E p = 610 keV resonance could be obtained due 
to the obscuring peak at E 1 — 6129 keV from 19 F con- 
tamination. In order to estimate the contribution from 
the 22 Na(p, 7) resonance, the spectral line shape was de- 
convolved into the contribution from 19 F and from 22 Na. 
This was done by comparing the on-resonance line shape 
with the sum of the off-resonance line shape and a shape 
representing the 6112-keV gamma ray. The latter shape 
was estimated from a normalized and shifted peak at 



F. Verification of Experimental Method 

In order to verify our technique, targets of 23 Na were 
implanted using the ion source on the injector deck at the 
low-energy end of the University of Washington accelera- 
tor, and known 23 Na(p, 7) resonances were measured un- 
der the same conditions as the 22 Na(p, 7) measurements. 
A 20-nm layer of chromium was also evaporated on the 
surface, similar to the main 22 Na targets. 

The 23 Na(p,7) resonance at E p = 512 keV has a re- 
ported strength of 91.3 ± 12.5 meV and is the recom- 
mended reference resonance for this reaction [49]. We 
measured this with a target containing (5.8 ± 0.9) xlO 15 
atoms, as determined by implantation charge integration 
with a measured correction for the sputtering of positive 
ions. The excitation function is shown in Fig. 1161 Using 
the 10810-keV gamma ray with a branch of 71% [50(, we 
determined its total resonance strength to be 79 ± 17 
meV, after applying the correction described below. 

Inspection of the raster plot indicated a small part of 
the target was missed when using our standard raster. 
To determine the missing fraction, the target was illumi- 
nated with a larger raster amplitude that fully encom- 
passed the target atoms. Comparing the full target area 
to the fractional area covered with the standard raster, 
the amount missed was estimated to be ~ 12%. This cor- 
rection was applied to the resonance strength. This was 
not necessary with our 22 Na targets, as their respective 
raster plots indicated the beam covered the entire active 
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area. To estimate a range for the effect of the observed 
23 Na target misalignment, simulations were performed, 
as described in Sec. Ill F 2| with an offset of 1.5 to 2.0 
mm and a variety of target distributions. The uncer- 
tainties associated with these calculations dominate the 
overall uncertainty quoted above. 



IV. RESULTS AND DISCUSSION 

Results are shown in Tables HI1 Hill and IIVI The reso- 
nance energies are summarized in Table UH the gamma- 
ray branches and partial strengths for each resonance 
are summarized in Table IIII1 and the final resonance 
strengths are summarized in Table IPV1 

The resonance energies determined from each method 
are shown in TableHIl and agreement between methods is 
quite good. The adopted energy is the weighted average 
of the two results. We find energies that agree with pre- 
viously reported values [HHIij], and we have improved 
the uncertainties on the energies. 

Table UlTl shows our excitation energies, gamma-ray en- 
ergies, branches, and partial strengths. Also included is a 
comparison with previous branches. Our branches are in 
agreement with the previous direct measurements [24l . [25[ 
for E p = 288 and 610 keV. For the 6112-keV branch of the 
610-keV resonance that could not be resolved from the 
19 F contamination, we found the upper limit to be 28%, 
which is consistent with the value of 20.0 ± 1.8% mea- 
sured directly by Stegmuller et al. [25[ . We have adopted 
the branch of Ref [25j to extract its partial strength be- 
cause it is consistent with ours and is more precise. For 
the previously established branches of the 454-keV reso- 
nance, our branches agree with Jenkins et al. [26j j , and we 
have improved upon their uncertainties. A new branch 
has also been identified with E 1 — 7566 keV. Although 
branches for the 213-keV resonance are in agreement with 
the earlier measurement [25j at the 2a level, an additional 
branch of ~ 11% has since been identified, and we agree 
with its currently established value f22j . 

For the observed resonances at E p = 213, 288, 454, 
and 610 keV, contributions to the total strength have 
been investigated for primary transitions to all levels up 
to the 6th excited state, also shown in Table UTTl For un- 
observed branches, limits were obtained via the method 
described in Sec. IIII El The branches with a final spin of 
1/2 are highly unlikely to be detectable based on angu- 
lar momentum considerations, so we attribute our non- 
zero result for E 1 = 5497 keV from the 288-keV reso- 
nance to statistical fluctuations. Upper limits on partial 
strengths from spin 1/2 levels have not been included 
in the strength uncertainties or in the branches. Other 
branches where only an upper limit could be determined 
were used to directly increase the upper uncertainty on 
the total strength but do not affect the central value of 
the total resonance strength. 

All finite total strengths are significantly larger than 
those previously reported [24| HH , as shown in Table IIVI 



The 213-keV resonance is stronger by a factor of 3.2, 
which includes a factor of 2.8 that we observe with re- 
spect to the same decay channels observed in Ref. [25j . 

Iacob et al. [22| attempted to extract a resonance 
strength for E p = 213 keV by compiling data from several 
different sources. Using the beta-delayed proton branch 
from Tighe et al. [20| , the proton-to-gamma-ray branch- 
ing of Perajarvi et al. [2l|, and the lifetime of Jenkins et 
al. [Ill, Iacob et al. claimed a total strength of 2.6 ± 0.9 
meV, in contrast to our value of 5.7±q;| meV. We con- 
sider the value of Iacob et al. to be unreliable. It is based 
on the observation of a /3-delayed proton peak very close 
to detection threshold that would have been difficult to 
disentangle from noise. Indeed, Saastamoinen et al. [5l| 
have shown that the /3-delayed proton intensity deduced 
by Tighe et al. [20| was too large by a significant amount. 
In addition, Iacob et al.'s value is based on the lifetime of 
the state, for which further documentation has not been 
published. 

We have set upper limits on the 198- and 232-keV po- 
tential resonances. For E p = 198 keV, based on a search 
for the 5055-keV branch, we observe a possible presence 
in our fits with a non-zero strength at slightly over la. 
Jenkins et al. [26] suggested that this resonance could 
have a strength as high as 4 meV, whereas our 68% C.L. 
upper limit is a factor of 8 smaller. The resonance at 
E p = 232 keV also was not observed, based on our search 
for the two main branches reported by Ref. [26| . This res- 
onance has been demonstrated to be the isospin analogue 
of the T = 3/2 23 Al ground state [HHH so its width for 
proton decay to T — 22 Na is significantly suppressed 
by isospin conservation, and thus it is expected to have 
a relatively small (p, 7) resonance strength. Our direct 
measurement is consistent with this expectation. 

Jenkins et al. 26] proposed a new level that corre- 
sponds to E p = 209.4(17) keV, which produces gamma 
rays with energies of 5729.1(11) and 5067.1(11) keV 
and branches of 33(6)% and 66(8)%, respectively. The 
gamma ray at 5067 keV is very close in energy to the 
5055-keV gamma ray from the possible resonance at 198 
keV, as they produce the same final state. We have inves- 
tigated this resonance, but due to the size of the energy 
windows necessary on each peak (~ 20 to 45 keV), it 
is unclear from which potential resonance these possible 
gamma rays originate. Contrary to Jenkins et al., Iacob 
et al. |22j attributed the gamma ray at 5729 keV to the 
resonance at 213 keV. Our branching ratios in Table [TTTl 
are consistent with Iacob et al. and not with Jenkins et 
al. An upper limit of 0.40 meV can be placed on the 
strength of the 209-keV resonance at the 68% confidence 
level, and it is possible that part of the contribution orig- 
inates from the potential resonance at 198 keV. 

Seuthc et al. 24] also measured resonances at E p = 
503, 740, and 796 keV that we did not investigate. These 
resonances do not play an important role in the nova 
scenarios we consider here, but we nevertheless included 
our estimate for their contributions in our calculation of 
the thermonuclear reaction rate. We assumed that the 
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TABLE II: Resonance energies. Values are extracted from the lab proton energy, E p , via the excitation function, correcting 
for energy loss E\ OBB , and from the gamma-ray energy, E 1 . Statistical uncertainties are given for each individual branch, and 
combined branches include systematic uncertainties. The adopted value is the weighted average of both methods. All values 
are in keV. 



JT i e V 1U Lla 


Present 


Hip 


E^ 


E p (from excitation function): 
-Bioss branch combined branches 


E p (from E-y): 
branch combined branches 


Adopted E p 


213.3 ± 2.7 


7332.7 ± 1.2 


5.1 ± 1.0 213.1 ± 3.0 213.1 ± 3.4 


213.6 ± 1.6 213.6 ± 1.6 


213.5 ± 1.4 


287.9 ± 2.1 


5140.6 ± 1.0 
5803.2 ± 1.3 


4.6 ± 0.9 286.5 ± 2.1 286.3 ± 2.3 
285.6 ± 4.2 


288.4 ± 1.7 288.7 ± 1.3 
289.1 ± 2.0 


288.1 ± 1.1 


457 ± 2" 


5300.1 ± 0.8 
5962.7 ± 0.8 


3.7 ± 0.8 452.8 ± 0.8 452.8 ± 1.1 
452.8 ± 1.1 


455.6 ± 1.7 455.7 ± 1.1 
455.9 ± 1.6 


454.2 ± 0.8 


611.3 ± 1.8 


8162.3 ± 0.9 
7711.2 ± 1.1 


3.2 ± 0.7 609.6 ± 0.9 609.0 ± 1.1 
608.1 ± 1.1 


611.0 ± 1.5 610.8 ± 1.2 
610.3 ± 2.0 


609.8 ± 0.8 



"Prom Ref. |25|| , unless otherwise noted. 
'From Ref. Q. 



relative strength of the resonances to the 454- and 610- 
keV resonances is correctly given by the result of Scuthe 
et al. Because we observe an average factor of 2.5 ± 0.5 
larger for the strengths of resonances observed here with 
respect to those in Ref. [Hf , we scaled the oj^/s in Ref. [HJ 
by this factor. 

It is surprising that our strengths are several times 
higher than those from the previously reported di- 
rect 22 Na(p, 7) measurements 0, HE|- However, those 
strengths were determined relative to resonance strengths 
from one experiment performed in 1990 24], also with 
implanted targets but with 2x the implantation energy 
and into a different substrate material. We rastered the 
beam over the entire target, whereas the measurement 
of Ref. [24| used a centered beam, which can degrade 
the target locally in ways that are difficult to charac- 
terize. The gamma-ray energy window in that experi- 
ment was several MeV for most of the data, although 
two resonances were measured at peak yield with a high 
resolution detector. Our energy windows were always 
narrow, with a maximum of a few tens of keV to incor- 
porate only the relevant peak. Along with a cosmic-ray 
anticoincidence system, our method employed full exci- 
tation functions integrated over proton energy and was 
largely independent of stopping-power estimations. Our 
method required knowledge of only the total number of 
target atoms, determined from the target activity, and 
the requirement that the beam covered all target atoms, 
which we could monitor with the information from the 
raster. The price paid for eliminating the dependence 
on the target distribution was the difficulty of determin- 
ing the beam density. This was accomplished experi- 
mentally using a 27 A1 coin target, and systematic effects 
were carefully considered. In addition, our detector effi- 
ciency was determined with two radioactive sources and 
two resonances in the 27 Al(p, 7) reaction on two different 
substrates, spanning a gamma-ray energy scale of 1.3 to 
11 MeV, in addition to calculation with detailed simula- 
tions. Furthermore, the validity of our method has been 
confirmed with the 23 Na(p, 7) measurement. A target 



was implanted, and result for the resonance strength are 
within 89 ± 24% of the currently accepted value. Al- 
though the error on this quantity is not negligible, it is 
not a factor of ~ 3 we observe in the 22 Na(p, 7) resonance 
strengths. Thus we are confident in the absolute values 
we have obtained for the strengths. 

V. ASTROPHYSICAL IMPLICATIONS 

Using the resonance energies and total strengths shown 
in Table IIV1 we calculated the contributions to the 
22 Na(p, 7) 23 Mg thermonuclear reaction rate under the 
narrow-resonance formalism. Monte Carlo methods were 
employed, assuming symmetric (asymmetric) gaussian 
distributions for all measured resonances with symmetric 
(asymmetric) uncertainties. For the proposed resonances 
at E p = 198 and 232 keV, their distributions were taken 
to be the curves shown in Fig. Q3] i? CC ntrai is defined as 
the 50% quantile of the distribution of rates at a given 
temperature, and i? U pper and i?i OW cr are the 16% and 84% 
quantiles, respectively. 

The proposed resonance at E p — 209 keV has not 
been included in the reaction rate, as its previously es- 
timated strength [2(| is so weak that its contribution 
should be negligible. In addition, the upper limit this 
work sets is conservative because of potential contribu- 
tions from nearby resonances to the excitation function. 
The strengths of the resonances at E p = 45 and 70 keV 
are too low to be measured directly due to the Coulomb 
barrier but have been included in the calculation of the 
reaction rate derived from R ef. [28l| , based on the ( 3 He,<f) 
spectroscopic factors of Ref fl9j . 

In Fig. [T7] (a), we show the individual contributions 
to the thermonuclear rate for 22 Na(p, 7), and panel (b) 
shows the relative contributions of selected resonances. 
As one can see, any contribution from the resonance at 
E p = 198 keV is less than that from the resonance at 
E p = 213 keV for all temperatures of interest to novae. 
Therefore, the resonance at E p = 198 keV does not dom- 
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TABLE III: Branches and partial strengths of 22 Na(p,7). E x and Et are the excitation and transition energies, respectively. 
The value given for the branch was determined from the sum of all data for a particular resonance. Systematic errors are 
not included in the partial strengths. 68% confidence levels are given for all data, and both an upper limit and measurement 
have been given for E v = 198 keV. Total strengths are shown in Table ITVl All upper limits for partial strengths were derived 
assuming their excitation function was the same shape as for the strongest branch. 

















Branches (%)" 






TP b 

E x 


E p 


IP b 

E t 


E f 






Previous 




Present 


Impartial 


(keV) 


(keV) 


(keV) 


(keV) 


If 


Ref. [24,1 


Ref. [25] d 


Ref. [26] 


Ref. [22] d 




(meV) 


7770. 2±1. 4 


198 


5055 


2715 


9/2 + , 5/2+ 




- 


58±8 


- 


* 


< 0.30 (0.20+0^3) 






2317 


5453 






- 


42±7 


- 






7782. 2±1. 2 


209 


5730 


2052 


7/2+ 


- 




33±6 


- 


- 


T 

J 






5067 


2715 


9/2+, 5/2+ 


- 


- 


66±8 


- 


- 


< 0.26 


7784. 7±1. 2 


213 


7785 





3/2+ 


- 


< 29 


_ 


3.1±2.0 


- 


< 0.09 






7334 


451 


5/2+ 


- 


100 


100 


80.8±3.6 


89.4±5.3 


5.1 ±0.5 






5732 


2052 


7/2+ 


- 


< 29 


- 


16.2±3.4 


10.6±5.3 9 


0.6±0.3 9 






5426 
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1/2+ 


- 


< 29 


- 


- 


- 
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- 


- 


- 
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- 


- 


- 
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- 


- 
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- 
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288 
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< 4.3 


- 


- 


- 
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5/2+ 
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- 
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- 
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- 


- 
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- 
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- 


- 


- 


- 


- 
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- 
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- 


- 
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- l 
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- 






- 
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5.3±1.7 






5963 
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7/2+ 






29±12 




43.6±1.2 


74.7±5.2 






5656 


2359 


1/2+ 












< 1.1 






5301 


2715 


9/2+, 5/2+ 


100 




71±16 




51.9±1.2 


85.9±5.5 






5244 


2771 


1/2" 












< 0.8 






5107 


2908 


(3/2,5/2)+ 












J 


8163.9±0.8 


610 


8164 





3/2+ 


65±5 


65.0±2.3 






61.3±1.8 


376±16 






7713 


451 


5/2+ 


19±2 


15.0±1.8 






18.6 ±1.3 


97±16 






6112 


2052 


7/2+ 


16±2 


20. Oil. 8 


100 




(20.0±1.8) fc 


(118±13)' 






5805 


2359 


1/2+ 




< 1-4 








< 24 






5449 


2715 


9/2+,5/2+ 




< 1-4 








< 11 






5393 


2771 


1/2- 












< 18 






5256 


2908 


(3/2,5/2)+ 












< 18 



a We assume the sum of all observed branches adds up to 100%. 
'Derived from our results in TablelTll Otherwise from NNDC [45j . 
c Branches from upper limits are calculated as the partial strength 
relative to the total observed strength. 
^Converted finite values into percents. 

^Partial strengths cannot be determined. See Table ITVl for upper 
limits on total strengths. 

■^Any contribution from this branch has been attributed to E p = 
213 keV. 

'Imposed the restriction that the shape of the excitation function 
must be the same as for the strongest branch. 

^Because If = 1/2, this state is highly unlikely to have a de- 
tectable value so we attribute this value to statistical fluctuations 
and do not include it in i^7 tota i. 

'Could not be determined due to a small percentage of pulser 
counts sorted into this energy window. 

3 Could not be determined due to 19 F background. 

*The value from Ref. [2511 is used, as this gamma ray was obscured 
by 19 F background in our measurements. 

'Estimated from branch. 
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TABLE IV: 22 Na(p,7) resonance strengths. Upper limits are 
at the 68% confidence level. Both an upper limit and measure- 
ment have been given for E p = 198 keV. The total strength 
was equal to the sum of partial strengths given in Table ITTT1 for 
all resonances, with the exception of E p — 198, 209, and 232 
keV where branches from Refs. [22l.[26| must be used. If only 
an upper limit has been given to a particular branch in Ta- 
ble lHIl then its value directly increased our upper uncertainty 
only, not the central value. 





Present 


Previous 


Present 


rp lab 


£p m (keV) 


k>7total (meV)° 


W7total (meV) 


45 


43.1 ± 1.7" 


(7.1 ± 2.9)xl0- 146 


(7.1 ± 2.9)xl0" i4 '' 


70 


66.6 ± 3.0 6 


(5.1 ± 2.1)xl0" loi> 


(5.1 ± 2.1)xl0" 106 


198 


189.5 ± 1.8 C 


< 4 C 


< 0.51 (0.34t°;i ) 


209 


200.2 ± 1.6 C 


(5 x 10~ 2 ) c 


< 0A0 d 


213 


204.1 ± 1.4 


1.8 ± 0.7 e 


c 7 +1.6 


232 


221.4 ± 2.3 C 


2.2 ± 1.0 C 


< 0.67 


288 


275.4 ± 1.1 


15.8 ± 3.4 


39 ± 8 


454 


434.3 ± 0.8 


68 ± 20 


166 ± 22 


503 


481 ± 2° 


37 ± 12 


93 ± 36 / 


610 


583.1 ± 0.8 


235 ± 33 


59lt^ 39 


740 


708 ± 2" 


364 ± 60 


913 ± 174 ; 


796 


761 ± 2° 


95 ± 30 


238 ± 79 / 



"From Ref. [2^, |25[ . 

'From or derived from Ref. [2g| . 

c From Ref. 

d Not included in the present reaction rate. 

E Actual value measured was 1.4 meV, inflated to 1.8 meV to ac- 
count for possible unknown branches. 
^Scaled from Ref 24]. Not measured in this work. 
9 Includes estimated partial branch from = 6112 keV. 
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FIG. 17: (color online) Thermonuclear 22 Na(p, 7) 23 Mg re- 
action rate as a function of temperature. Panel (a) shows 
contributions to the reaction rate from individual resonances 
labeled by E v (in keV), based on present measurements. 
Hatched regions represent 68% confidence levels (although the 
uncertainties are so small they may appear as solid lines), and 
dashed lines with arrows are 68% confidence level upper lim- 
its. Panel (b) shows the fractional contributions of selected 
resonances to the total rate, as calculated using the resonances 
illustrated in panel (a). 



inate the reaction rate in this region, as Jenkins et al. [261 ] 
proposed it might. Rather, the resonance at E p = 213 
keV makes the most important contribution. Also, at the 
higher temperatures around 0.4 GK, the contribution of 
the 288-keV resonance becomes significant. The other 
resonance contributions are effectively negligible at nova 
temperatures. 

Fig. [18] illustrates the total reaction rate relative to 
previous direct measurements [24| [H| , showing that our 
rate is inconsistent with previous work at all tempera- 
tures. The dashed line and hatched region in this figure 
do not include the distributions for proposed resonances 
at E p — 198 and 232 keV. The 232-keV resonance makes 
an insignificant contribution at all temperatures of inter- 
est to novae. However, although the 198-keV resonance 
was unobserved and an upper limit has been placed on 
its strength, including its potential contributions to the 
reaction rate has a non-negligible effect at low temper- 
atures. Including the proposed resonances increases the 
upper limit of the reaction rate to the dot-dashed line 
shown in Fig. [18] Because of this difference, we have cal- 
culated the reaction rate for each of these two separate 
cases, and the values are shown in Table IVl 

We can anticipate the general ramifications of the new 
rate on expected nucleosynthesis of 22 Na in ONe novae 
using post-processing network calculations because the 
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Na(p,7) Mg reaction rate as a function of temperature. 
Hatched areas represent 68% C.L. error bands and dashed 
lines represent the central value, each relative to the previous 
central value. Both previous and present resonance strengths 
are listed in Table ITVl Including the distributions for the res- 
onances at E p = 198 and 232 keV increases the present upper 
limit to the dot-dashed line. The upper limit for the 198-keV 
resonance from Jenkins et al. [26j| has not been included in 
the previous rate. 
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TABLE V: Table of the thermonuclear 22 Na(p, 7) 23 Mg reaction rate, R, as a function of temperature, T, determined using 
the energies and strengths given in Table ITVl Two distributions of rates were calculated, one of which includes the probability 
density functions from resonances at E p = 198 and 232 keV and one which does not. i? ccn t r ai is the 50% quantile of the 
distribution of rates at a given temperature, and -Ru PP cr and i?i OWGr are the 16% and 84% quantiles, respectively. The units of 
the rate are cm 3 mol _1 s _1 . 
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8.1xl0+ 02 


7.6xl0+ 02 


6.9xl0+° 2 


8.3xl0+ 02 


1.5 


2.2xl0+° 3 


2.1xl0+ 03 


2.4xl0+ 03 


2.2xl0+ 03 


2.1xl0+° 3 


2.4xl0+ 03 


2.0 


3.9xl0+° 3 


3.6xl0+ 03 


4.2xl0+° 3 


3.9xl0+ 03 


3.6xl0+° 3 


4.2xl0+ 03 



a "198, 232" denotes resonances at E p = 198 and 232 keV. 



TABLE VI: 22 Na yields from novae on ONe white dwarfs of 
various masses computed with SHIVA 0]. The total mass 
ejected is given for each case, along with the mass fractions of 
22 Na obtained using the central value of the previous [241. 125I] 
and both present 22 Na(p, 7) 23 Mg reaction rates compiled in 
Table FVl The factor given is the ratio of previous to respective 
present amounts of 22 Na ejected. 





1.15M 


1.25M© 


1.35M a 


Mojcct (g) 


4.9xl0^ s 


3.8xl0 2S 


9.0xl0 2 ' 


Previous 6 


1.6xl0 -4 


1.9xl0~ 4 


5.9xl0~ 4 


Present c 


8.8xW~" 


l.lxlO" 4 


4.1xl0" 4 


Factor 


1.8 


1.8 


1.4 


Present d 


7.8xl0~ 5 


l.OxlO" 4 


3.8xl0" 4 


Factor 


2.0 


1.9 


1.5 



"Statistically, there should be more novae of 1.15Mq and 1.25Mq 
than those hosting 1.35M0, due to the stellar mass function of the 
progenitors. 

'Rate derived using energies and strengths from Ref. [24l. |25||. 
c Not including "198, 232" from Table |V] 
including "198, 232" from Table [V] 



total energy generation is not affected appreciably by the 
22 Na(p, 7) reaction. Based on the one-zone calculations 
of Ref. a specific model indicates that the produc- 



tion of 22 Na is related inversely to the 22 Na(p, 7) reaction 
rate. The authors of Ref. [49| also varied the rate by a 
factor of two using various one-zone models of novae to 
extract the effect. Using the information given in these 
references and our change in the reaction rate, the es- 
timated abundance of 22 Na in novae is expected to be 
reduced by factors of 2 to 3 of what was previously ex- 
pected, depending on white dwarf mass and composition. 
This will directly affect the expected flux of 22 Na gamma 
rays observed using orbiting gamma-ray telescopes. 

In addition, the impact of our 22 Na(p, 7) reaction rate 
on the amount of 22 Na ejected during nova outbursts 
has been tested through a scries of hydrodynamic simu- 
lations: three evolutionary sequences of nova outbursts 
hosting ONe white dwarfs of 1.15, 1.25 and 1.35A/ Q 
have been computed with the spherically symmetric, La- 
grangian, hydrodynamic code SHIVA, extensively used in 
the modeling of such explosions (see Ref. Q , for details) . 
Results have been compared with those obtained in three 
additional hydrodynamic simulations, for the same white 
dwarf masses described above and same input physics 
except for the 22 Na(p,7) rate, which was derived from 
Refs. 0, HH . The network used for additional reaction 
rates is the relevant subset of that used in Ref. [52j . The 
estimated 22 Na yields (mass- averaged mass fractions in 
the overall ejected shells) are listed in Table I VII which 
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clearly shows that the impact of the central value of new 
rate roughly translates into lower 22 Na abundances by 
a factor up to ~ 2 with respect to previous estimates. 
This, in turn, directly affects the chances to potentially 
detect the 1275-keV gamma-ray line from 22 Na decay, de- 
creasing the maximum detectability distances by a factor 
~ 1.4. The inclusion of the 198- and 232-keV distribu- 
tions in the rate does not appreciably alter this factor. 

The results from one-zone post-processing network 
calculations and full hydrodynamic simulations using 
SHIVA are complementary. The post-processing ap- 
proach mimics the processes that occur in the deepest 
envelope layers, whereas the hydrodynamic simulations 
average the yields over all ejected shells. Convection also 
plays a critical role, supplying freshly, unburned mate- 
rial from external shells into the innermost one (and vice 
versa), and these effects cannot be simulated in a post- 
processing framework. As a result of the more realistic 
physics in the hydrodynamic model, the composition of 
the innermost shell is diluted by the compositions of the 
outermost ones. On the other hand, the post-processing 
calculations cover various nova models, a wider range of 
nova masses and compositions, and show that the corre- 
lation between the 22 Na(p, 7) rate and 22 Na production is 
robust even when other reaction rates are simultaneously 
varied. It seems reasonable to assume that the magnitude 
of the dilution from the hydrodynamic models (m 25%) 
applies generally to all of the post-processing results. 

VI. CONCLUSIONS 

We have measured the resonance strengths, energies, 
and branches of the 22 Na(j?, 7) 23 Mg reaction directly and 
absolutely. Our method improved upon past measure- 
ments in several ways. The use of integrated yields 
makes the results independent of absolute stopping power 
calculations and is far more robust than using peak 
yields. We also utilized isotopically-pure, implanted tar- 
gets that demonstrated nearly zero loss during bombard- 
ment. HPGe detectors exhibit excellent energy resolu- 



tion, providing the ability to use narrow energy win- 
dows, and anticoincidence shields enabled suppression of 
the cosmic ray background. Absolute detector efficiency 
was also vital, which we determined by fusing measure- 
ment and simulation. Finally, the rastering of the beam 
across the target not only aided in maintaining target 
integrity, but also removed the requirement of detailed 
knowledge of the target distribution. A determination 
of the beam density was mandatory and was ascertained 
by both measurement and modeling. As a consequence 
of the aforementioned points, our results should be sub- 
stantially more reliable than previous measurements. 

By exploiting these advantages, our measurement has 
shown that four previously measured resonance stren gths 
are 2.4 to 3.2 times higher than previously reported [24|, 
I25I ]. Jenkins et al. also proposed that a new 22 Na(p, 7) 
resonance with E p — 198 keV could dominate this reac- 
tion rate and have consequences for novae [26]. We have 
demonstrated that this is not the case, and that the main 
contributions arise from the resonance at E p = 213 keV. 
As a result of the higher resonance strengths, the esti- 
mated flux of 22 Na gamma rays from novae is expected 
to be about a factor of 2 less than what was previously 
expected, determined by using both post-processing net- 
work calculations and hydrodynamic simulations. The 
lack of observational evidence of 22 Na in the cosmos is 
consistent with the previous reaction rate; however, the 
present rate makes detection ~ 1.4 times more difficult 
to detect. 
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